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摘 =: 利用 2015 年 8 月 一 2016 年 7 月 在 印度 河上 游 流 域 Bagrot 山谷 降水 稳定 同位 素 (3 "0 和 
8D) 观 测 结果 以 及 当地 气象 资料 ,利用 同位 素 示 踪 及 统计 分 析 方 法 ,并 结合 HYSPLIT 模型 ,对 研究 
区 降水 稳定 同位 素 变化 特征 大 气 水 线 以 及 水 汽 来 源 进 行 了 分 析 。 结 果 表明 :观测 期 间 Bagrot 山谷 
降水 稳定 同位 素 的 季节 变化 明显 ,8”0 与 8D 秋冬 季 偏 低 , 春 夏季 偏 高 , 且 与 气温 变化 一 致 ,存在 显 
著 的 温度 效应 ,而 降水 量 效 应 不 明显 。 而 上 且 发 现 研究 区 局 地 大 气 水 线 截 距 和 和 斜率 均 低 于 全 球 的 , 反 
映 了 降水 过 程 中 云 下 二 次 茹 发 作用 较为 强烈 ,因此 ,不 同 的 降水 形态 导致 该 研究 区 局 地 大 气 水 线 的 
斜率 和 截 距 不 同 。 当 液态 降水 (降雨 ) 发 生 时 ,由 于 在 较为 干旱 的 气候 环境 下 ,雨滴 在 降落 的 过 程 
中 受到 二 次 蒸发 相对 较 强 , 使 得 局 地 大 气 水 线 的 斜率 和 截 距 偏 低 ; 而 当 固 态 降水 ( 降雪) 发生 时 ,由 
Ch 于 温度 较 低 , 受 再 循环 水 汽 和 二 次 燕 发 的 影响 较 小 ,导致 局 地 大 气 水 线 的 斜率 和 截 距 均 偏 商 。 
Tr Bagrot 山谷 及 其 周边 地 区 ,从 南 到 北 局 地 大 气 水 线 的 斜率 相差 不 大 ,而 其 截 距 总 体 上 随 着 纬度 升 高 
而 降低 ,可 能 与 云 下 二 次 茹 发 导致 稳定 同位 素 发 生 的 不 平衡 分 馏 逐 渐 强烈 有 关 。 通 过 Bagrot 山谷 
上 点 降水 稳定 同位 素 观测 结果 并 结合 HYSPLIT 模型 的 后 向 追踪 ,研究 还 发 现 ,研究 区 全 年 主要 受 
= 西风 环流 以 及 局 地 环流 的 影响 。 但 与 研究 区 以 北 的 临近 站 点 ( 慕 土 塔 格 、 和 田 等 ) 相 比 有 所 不 同 ， 
a 由 于 Bagrot 山谷 位 置 更 靠 南 , 其 仍然 偶尔 受到 来 自 南方 的 海洋 性 水 汽 影响 。 这 一 研究 结果 可 能 对 
¢ 该 地 区 树 轮 稳定 同位 素 记 录 的 解 译 具有 一 定 的 指示 意义 。 

r= 关键 词 : 印度 河上 游 ; 降水 稳定 同位 素 ; 大 气 水 线 ; 水 汽 来 源 
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被 称 为 水 的 “指纹 ”的 稳定 同位 素 (8D 和 流域 水 循环 过 程 以 及 区 域 气候 变化 的 理解 ,同时 也 
8^0) ,能 非常 敏感 地 响应 气候 环境 的 变化 ,并 记载 GS ` 石 筹 , 树 轮 中 稳定 同位 素 记录 重建 古 气候 提 
着 水 循环 的 整个 历史 信息 。 大 量 研究 表明 ,降水 。” 供 重要 的 理论 基础 。 
稳定 同位 素 变化 主要 受 温度 、 降 水 量 、 相 对 湿度 、 气 Hie 1964 年 ,DANSGAARD fè h T ESR (d- 


压 等 不 同 气象 要 素 的 影响 ” ,在 不 同 海拔 .下 热 面 
以 及 莹 发 与 凝结 程度 的 不 同 也 会 导致 降水 稳定 同位 


excess 或 简写 为 4) 的 概念 ,将 其 定义 为 d= 5D - 
88°07), d 反映 了 水 体 在 动力 分 馏 过 程 中 稳定 同 


素 发 生 改 变 ””。 另 外 ,降水 稳定 同位 素 还 与 不 同 
的 水 汽 来 源 、 水 汽 传输 过 程 及 输送 方式 联系 密 
GS 。 因 此 ,降水 稳定 同位 素 研 究 有 助 于 加 深 对 
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位 素 变化 的 差异 ,主要 受 水 汽 源 地 的 温度 以 及 相对 
湿度 影响 ,同时 也 与 降水 过 程 中 的 气象 条 件 有 关 "" 。 
已 有 研究 表明 ,全球 降 水 d 平均 值 约 为 10%o” ,而 相 


作者 简介 : 


王 邮 凡 〈1991 - ) , 男 ,河南 安阳 人 ,在 读 硕士 研究 生 , 主 3 


通讯 作者 : 张 寅 生 . E-mail; yszhang@ itpcas. ac. cn 


研究 方向 为 地 表 水 文 过 程 . E-mail; wangyefan@ itpcas. ac. cn 


对 于 10%o 较 高 的 d 值 代表 水 汽 源 地 气候 较为 干燥 ， 
水 汽 燕 发 过 程 中 受到 了 强烈 的 动力 分 馏 的 作用 , 较 
低 的 a 值 则 代表 了 较为 湿润 的 海洋 性 气 团 带 来 的 水 
汽 "1。 另 外 ,在 降水 过 程 中 , 云 下 二 次 藻 发 会 加 剧 
同位 素 动力 分 馏 作用 ,也 会 导致 降水 中 d 值 明 显 降 
低下 。 因 此 ,利用 4d 可 以 为 研究 水 汽 源 地 、 水 循 
环 以 及 局 地 气候 条 件 提供 了 一 定 的 依据 ”|。 

全 球 大 范围 的 降水 稳定 同位 素 观测 活动 可 以 追 
WEJ 1961 年 ,国际 原子 能 机 构 (IAEA) 与 世界 气象 
组 织 (WMO) 合 作 在 全 球 范 围 内 建立 了 降水 稳定 同 
位 素 观 测 网 络 (Global Network of Isotopes in Precipi- 
tation, GNIP) ,先后 在 全 球 120 多 个 国家 和 地 区 建 
立 了 1 000 多 个 降水 观测 站 点 ,以 期 对 全 球 降水 中 
稳定 同位 素 进行 长 期 监测 。 中 国 对 降水 同位 素 
的 研究 虽然 开始 较 晚 ,但 发 展 较为 迅速 ,无 论 是 在 全 
Fa BE E? ,还 是 地 区 流域 尺度 上 -都 取得 
了 很 大 的 进展 。 尤 其 是 从 1991 年 开始 ,中国 科 学 院 
竺 青藏 高 原 上 启动 了 降水 稳定 同位 素 观 测 计划 
(Tibetan Plateau Network of Isotopes in Precipitation , 
TNIP) 77-9 SAA 40 多 个 站 点 。 以 这 些 站 
点 观测 资料 为 基础 ,YAO 等 揭示 了 青藏 高 原 降 水 稳 
定 同位 素 表现 出 季风 模 态 .过 渡 模 态 和 西风 模 态 等 
三 种 形态 "| ;其 中 在 主要 受 西风 影响 的 青藏 高 原 北 
部 ,降水 8^0 主要 表现 为 温度 效应 (3 "0 与 温度 呈 
正 相 关 关系 ) ;在 青藏 高 原 中 部 过 渡 区 ,在 长 时 间 尺 
度 上 ,降水 S" O 受 温度 变化 控制 ,表现 为 温度 效应 ， 
但 夏季 印度 季风 降水 会 削弱 其 与 温度 的 关系 ;而 在 
主要 受 印 度 季风 影响 的 青藏 高 原 南部 ,降水 3 0 FE 
要 表现 为 降水 量 效应 (3 O 与 降水 量 呈 反 相 关 关 
系 ) “1。 
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了 部 分 成 果 , 但 主要 集中 于 基于 单 次 考察 ,缺乏 连续 
性 , 故 对 水 体 稳定 同位 素 的 季节 变化 研究 还 不 够 深 
入 , 且 主要 侧重 于 地 表 水 体 稳 定 同 位 素 研究 “ =n 
地 表 水 对 于 捕捉 水 汽 来 源 温度 .降水 量变 化 以 及 蒜 
发 ,对流 .凝结 等 天 气 过程 不 够 敏感 并 且 具 有 兆 后 
性 。 相 对 而 言 , 降 水 稳定 同位 素 研究 能 够 有 效 地 弥 
补 这 一 不 足 。 连 续 而 系统 地 观测 降水 稳定 同位 素 ， 
对 印度 河上 游 水 汽 来 源 变化 及 其 水 循环 过 程 研 究 具 
有 积极 的 意义 。 在 此 背景 下 ,本 研究 在 印度 河上 游 
地 区 Bagrot 山谷 连续 一 年 开展 了 降水 稳定 同位 素 观 
测 ,并 结合 当地 气象 资料 ,探讨 了 研究 区 降水 稳定 同 
位 素 的 季节 变化 特征 ,并 对 其 水 汽 来 源 进 行 了 初步 
分 析 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

降水 样品 收集 站 点 Bagrot 山谷 (74. 56° E, 
36.03? N, 海拔 2 390 m) 位 于 巴基斯坦 北部 印度 
河上 游 (图 1) ,处 在 喀 喇 昆仑 山脉 南 坡 ,山谷 总 长 度 
为 32 km。 据 1994 一 2009 年 山谷 内 海拔 2 310 m 处 
气象 观测 资料 ,山谷 内 年 平均 温度 为 11.6 % ,年 平 
均 降 水 量 为 293.6 mm ,月 平均 降水 和 温度 (图 1b) , 
主要 分 布 在 春 夏 二 季 , 占 全 年 降水 的 70%。 山 谷 内 
由 于 降水 较 少 ,植被 稀 下 ,表现 为 典型 的 半 干 旱 气 候 
SED. 
1.2 样品 收集 与 分 析 

在 2015 年 8 月 一 2016 年 7 月 一 个 完整 水 文 年 
共 收 集 83 个 降水 样品 ,包括 71 个 液态 降水 样品 以 
及 12 个 固态 降水 样品 。 对 于 液态 降水 ,在 每 次 降水 


位 于 青藏 高 原 西 缘 的 印度 河上 游 , 地 处 西风 区 
和 季风 区 交界 处 ,西风 与 印度 季风 如 何 影 响 其 区 域 
内 水 循环 变化 广 受 关注 。KARIM 等 人 利用 河水 稳 
定 同位 素 计算 得 出 印度 河上 游 有 64% ~72% 水 汽 
源 于 地 中 海地 区 , 由 西风 环流 带 来 “1 。 类 似 的 ， 
SHARMA 指出 印度 河上 游 由 西风 环流 带 来 的 水 汽 
为 51%。 二 者 研究 结果 之 所 以 有 差别 是 因为 他 们 
研究 的 具体 区 域 不 同 :前 者 研究 区 域 位 于 印度 河上 
游 在 巴 控 克 什 米尔 地 区 ,而 后 者 主要 聚焦 于 印度 河 
上 游 在 印 控 克什米尔 地 区 。 相 对 而 言 ,西风 对 于 前 
者 的 影响 比 后 者 的 (更 靠 西南 ) 影 响 要 更 为 强烈 , 因 
此 由 前 者 的 研究 区 域 西风 环流 所 带 来 的 水 汽 含量 更 
大 号 ] 。 印 度 河 上 游 水 体 稳定 同位 素 研 究 虽 然 取 得 


事件 结束 后 ,立即 将 其 倒 入 15 mL 的 样品 瓶 密封 保 
存 ; 对 于 固态 降水 , 则 先 装 入 样品 袋 中 予以 封 紧 , 待 
其 在 室温 下 完全 融化 后 ,再 装 入 样品 瓶 中 密封 保存 。 
所 有 样品 在 进行 实验 室 测量 之 前 均 冷 冻 保 存 。 在 降 
水 样品 收集 同时 ,对 降水 事件 的 开始 与 结束 时 间 A 
温 、 降 水 量 进行 记录 。 

样品 的 分 析 工 作 在 中 国 科学 院 青藏 高 原 研 究 所 
环境 变化 与 地 表 过 程 重 点 实验 室内 用 Picarro- 
L2130i 液体 水 同位 素 光 谱 仪 上 进行 的 。 稳 定 同位 
素 的 比率 用 相对 于 维也纳 标准 平均 海洋 水 (V- 
SMOW2 ) 的 千 分 差 表 示 : 


8 (%o) = Coe 1) x1 000 (1) 
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图 1 印度 河上 游 位 置 及 降水 样品 收集 站 点 示意 图 (a. V X EEE URE a LN A ; 
b. Bagrot 山谷 多 年 温度 和 降水 月 变化 ) 


Fig.1 Location of Upper Indus Basin and the sampling site. (a. Sketch of the westerlies, Indian monsoon and 


Upper Indus Basin; b. Monthly variations in temperature and precipitation amount in Bagrot Valley) 


式 中 Res 和 Ry smow2 分 别 为 降水 样品 和 V-SMOW2 
的 稳定 同位 素 比 值 (D/ 孔 或 “0 和 0) ,仪器 测试 精度 
分 别 为 0.6%o (8D) 和 0. 1%o (3 0) 。 降 水 稳定 同 
位 素 的 日 月 和 季节 平均 值 为 降水 量 加 权 平 均 ,计算 
方法 如 下 : 


8^0,- Y Pia"0,/ ES (2) 


式 中 :8*0; 和 P; 为 次 降水 8^ 0. 以 及 对 应 的 降水 量 ， 
8^0, 为 日 月 降水 8^0 加 权 平 均值 。 

为 验证 稳定 同位 素 所 指示 的 降水 水 汽 来 源 ,本 
研究 选取 美国 海洋 及 大 气管 理 局 (NOAA ) 所 提供 的 
HYSPLIT 气 团 轨迹 模型 31 (http: // ready. arl. noaa. 
gov/HYSPLIT. php) ,利用 美国 国家 环境 预测 中 心 
(NCEP) 的 全 球 再 分 析 资 料 (ftp: // arlftp. arlhq. no- 
aa. gov/pub/archives/reanalysis ) ,以 降水 样品 收集 站 
点 为 示 踪 终点 ,采用 后 向 轨迹 法 ,追踪 观测 站 点 上 空 
1 500 m 气 团 120 h 内 的 后 向 传输 轨迹 。 


2 GRAD 


2.1 降水 稳定 同位 素 的 变化 及 其 与 气象 要 素 的 
关系 

图 2 和 表 1 展示 了 Bagrot 山谷 2015 年 8 月 一 
2016 4E 7 H d D5" O 温度 和 降水 量 的 日 变化 特 
fiE, 8^O 5 8D 变化 范围 分 别 为 - 22. 54%o ~ 
9.70%o 和 — 165. 69%o ~ 81.90%o, 年 加 权 平 均值 分 别 
为 -5. 92%o 和 — 40. 87%o。d 值 的 变化 范围 为 
-14.19%o ~21.52%o, 年 加 权 平 均值 为 6. 52%o, 低 
于 全 球 平均 值 10%o。 由 于 研究 区 降水 稳定 同位 素 
All d 值 季节 变化 较为 明显 ,本 研究 将 全 年 分 为 秋冬 
季 (9 月 ~ 翌年 2 月 ) 和 春 夏 季 (3 ~ 8 月 ) 进 行 分 析 。 

Bagrot 山谷 秋冬 季 与 夏秋 季 降 水 稳定 同位 素 和 
平均 值 差异 明显 ,秋冬 季 80 与 8D 加 权 平 均值 分 
别 为 -4.86%o 和 -33. 16%o, 春 夏季 8^ 0 5j 8D 加 
权 平 均值 分 别 为 - 10. 56%o 和 — 74. 48%o, 秋 冬季 
六 0 与 SD 明 显 低 于 春 夏季 的 。 具 体 表 现 为 ,秋季 
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表 1 6*O.6D 和 4 的 最 高 值 .最 低 值 以 及 加 权 平 均值 


Tab.1 The maximum-value, minimum-value and weight average value for O, 8D and d 
in different time scales ( whole year, autumn-winter, and spring-summer ) 
全 年 秋冬 季 春 夏 季 
8^0 8D / %o d 8^0 8D / %o d 8^0 8D / %o d 
最 低 值 -22.54 -165.69 -14.19 -22.54 -165.69 -7.36 -17.99 -144.14 -14.19 
最 高 值 9.70 81.90 21.52 7.50 52.64 18.04 9.70 81.90 21.52 
加 权 平 均值 -4.86 -33.16 6.52 -10.56 -74.48 9.96 -4.86 -33.16 5.72 
E 一 秋冬 季 一 一 上 一 一 者 夏季 在 后 续 研 究 中 结合 更 长 时 间 尺 度 的 降水 稳定 同位 素 
20 |, @ 过 量 气 观测 ,以 期 进行 深入 分 析 。9 月 23 日 温度 又 降 , 并 
| 且 对 应 的 降水 量 为 4. 8. mm, 水 汽 冷凝 较为 充分 , 降 
100 水 量 较 大 ,雨滴 在 降落 的 过 程 中 受到 的 二 次 蒸发 作 
0 x ”用 较 弱 ,因此 降水 8*0 的 值 也 相应 发 生 又 降 , 达到 
ToS 了 全 年 最 低 值 -22.54%o。 之 后 ,降水 880 与 8D f 
-200 。 。 随 着 温度 的 下 降 逐 渐 降低 。 
3 月 人 春 开 始 , 温 度 开 始 上 升 , 对 应 的 8^0 与 
8D 开始 波动 上 升 ,到 夏季 达到 最 高 值 。 从 图 3a ~ 
为 ”3c、 表 2 看 出 ,Bagrot 山谷 降水 8*0 与 8D 与 温度 之 
WES 间 均 存在 良好 的 线性 关系 ,在 全 年 ,秋冬 季 和 春 夏季 


0 
2015-08 2015-10 2015-12 2016-02 2016-04 2016-06 2016-08 
年 份 
图 2 2015 年 8 月 一 2016 年 7 月 Bagrot 山谷 
各 气象 要 素 和 同位 素 的 季节 变化 
Fig.2 Seasonal changes in meteorological elements, 


8D and 80 in Bagrot Valley from August, 2015 to July,2016 


8^0 与 8D 随 着 温度 的 下 降 开 始 波动 下 降 , 进 入 冬 
季 后 温度 达到 全 年 最 低 ,同时 对 应 的 8*0 与 5D 也 
达到 全 年 最 低 并 且 相 对 比较 稳定 ,降水 多 以 固态 形 
式 (降雪 ) 出 现 。 结 合 8s0 的 值 和 气象 数据 来 看 , 秋 
冬季 降水 较为 稀少 ,整个 秋冬 季 仅 有 21 d 发 生 降水 
事件 。 鉴 于 入 秋 以 后 9 月 份 的 降水 事件 ,9 月 20 日 
8^0 的 值 为 7.50%o, 远 高 于 其 前 后 的 9 月 7 日 以 及 
9 月 23 日 降水 80 的 值 ( -1.37%o 和 -22.54%o) 。 
人 研究 发 现 , 此 日 温度 为 15.3 % ,明显 高 于 9 月 7 日 
和 9 月 23 日 发 生 降水 事件 时 的 温度 (12.3 CC 和 
11.3 C) ,同时 对 应 的 降水 量 较为 稀少 ( 仅 为 1.1 
mm)。 因 此 ,9 月 20 日 降水 8^0 的 异常 高 值 , 可 能 
与 温度 升 高 以 及 细小 雨滴 在 降落 过 程 中 受到 明显 的 
二 次 莹 发 有 关 , 也 可 能 还 与 其 它 要 素 有 关 ,本 研究 将 


等 不 同时 期 ,二 者 均 呈 正 相 关 关 系 ,相关 系数 分 别 为 
0.73 .0.72 和 0.66 , 且 均 能 够 通过 已 < 0.01 的 显著 
性 检验 ,表明 Bagrot 山谷 降水 稳定 同位 素 存 在 显著 
的 “温度 效应 ”。 而 且 发 现 , 秋 冬季 降水 8^0 与 温度 
的 相关 系数 (0.72) 高 于 春 夏 季 的 (0. 66) ,这 说 明 
在 秋冬 季 温 度 对 于 降水 稳定 同位 素 有 着 绝对 的 控 
制 ,而 春 夏 季 温 度 对 降水 稳定 同位 素 的 影响 相对 减 
弱 。 但 Bagrot 山谷 降水 稳定 同位 素 与 降水 量 无 明显 
的 相关 关系 ,不 表现 出 “降水 量 效 应 ” (图 3d ~3e, 表 
2)。 这 与 前 人 的 人 研究 结果 , 即 内 陆地 区 降水 量 效 应 
不 明显 ,而 降水 8^ O 在 内 陆 中 高 纬度 地 区 更 易 受 到 
温度 的 影响 -是 一 致 的 。 在 青藏 高 原 西北 部 的 
zi uL; 518] 0^ EN EO) 、 和 田力 至 帕 米 尔 高 原 
的 塔什库尔干 RAE EO RA 
y 7*9 eh PEE Je 91 ,降水 稳定 同位 素 具 有 明 
显 的 温度 效应 ,主要 受 西 风 环流 以 及 局 地 水 循环 的 
影响 。 

Bagrot 山谷 d 值 全 年 季节 变化 表现 为 春 夏 低 ， 
秋冬 高 ,波动 范围 较 大 。 其 中 , 春 夏 季 d 值 的 变化 范 
围 为 : — 14. 19%o ~ 21. 52%o ,平均 值 为 5.72%o ,秋冬 
Æ d 值 的 变化 范围 为 : -7.36%o ~ 18. 04%o ,平均 值 
Ky 9. 96%o ,接近 全 球 的 平均 水 平 (10%o) 。 春 夏季 d 
值 小 于 10%o 的 次 数 (40 次 ) 远 远 多 于 d EKF 10960 
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图 3 不 同时 间 尺 度 下 Bagrot 山谷 降水 稳定 氧 氧 稳 定 同位 素 与 温度 和 降水 量 的 关系 


Fig.3 Relationship between O (8D) and temperature and precipitation in Bagrot Valley for 


#2 5°O 6D 与 温度 和 降水 量 的 线性 关系 
Tab.2 Correlation between stable isotopes (O and 8D) 
and meteorological factors (temperature and precipitation 
amount) at Bagrot Valley for different time scales 


( whole year, autumn-winter, and spring-summer ) 


80.T 8D -7 8/50-P 8D-P 
全 年 ”斜率 0.87 6.60 0.043 0.51 
截 距 -15.34 — -113.04 -6.37 -46.30 
r 0.73 0.71 0. 03 0. 05 
n 72 72 72 72 
P 0.00 0.00 0.80 0.70 
秋冬 季 斜率 1.04 8.04 0.201 1.963 
截 距 -17.36 — - 128.41 -11.82 - 86. 79 
r 0. 72 0. 72 0. 09 0. 11 
n 21 21 21 21 
P 0.00 0.00 0.71 0. 64 
ABE qp 0.74 5.63 -0.14 -0.85 
WIE -13.29 -98.24 -3.24 -23.36 
r 0.66 0. 62 0. 12 0. 09 
n 51 51 51 51 
P 0.00 0.00 0.39 0.52 


TE: 7 为 相关 系数 ,n 为 样品 数 ,P 为 显著 性 水 平 


的 次 数 (11 次 ) , 春 夏季 的 降水 量 占 到 了 全 年 降水 量 
的 81.3% , 春 夏 d 值 较 秋 冬季 低 ,并且 波动 幅度 较 
大 ,表明 该 时 段 内 水 汽 来 源 比较 为 复杂 。 尤 其 是 , 研 
究 发 现 , 春 夏 季 降 水 中 d 值 发 生 较 低 事 件 对 应 的 降 
水 量 都 很 小 ,雨滴 在 降落 的 过 程 中 经 历 了 明显 的 二 
次 蒸发 ,动力 分 馏 作用 较 强 ,从 而 导致 降水 中 a 值 出 
现 异常 低 值 。 秋 冬季 的 降水 量 只 占 全 年 降水 量 的 
18.7% , 而 秋冬 季 d 值 整体 维持 在 一 个 比较 高 的 值 
域 范围 内 ,这 可 能 与 降水 的 固态 形式 有 关 , 降 雪 事 件 
发 生 时 ,由 于 环境 温度 相对 较 低 ,因此 降雪 在 降落 过 
程 中 的 动力 分 馏 相 对 较 弱 ,从 而 使 得 其 d. 值 相对 较 
m " 
2.2 局 地 大 气 水 线 

降水 稳定 同位 素 中 8D 5j 8^0 之 间 的 线性 关系 
称 为 大 气 降水 线 。 在 1961 年 CRAIG 根据 来 自 全 球 
不 同 地 区 的 400 多 个 降水 .河水 及 湖水 样品 中 SD 
H 8"0 定义 了 全 球 大 气 水 线 (Global Meteoric Wa- 
ter Line ,简称 GMWL) ,这 一 规律 可 以 用 8D =85 0 
+10 来 表示 "1 。 随 后 ,ROZANSKI 对 全 球 206 个 站 
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点 长 期 观测 的 降水 8D 与 9^0 进行 分 析 , 进 一 步 修 
正 了 大 气 水 线 方程 , 即 8D = 8.178*O+10.35 Pl, 
局 地 大 气 水 线 (Local Water Meteoric Line ,简称 
LWML) 是 指 某 一 地 区 降水 中 8D 和 5"0 之 间 的 线 
性 关系 。 水 汽 从 源 地 蒸发 到 凝结 降水 的 这 一 过 程 
,影响 稳定 同位 素 分 馏 的 因子 有 所 不 同 ,因此 局 地 
大 气 水 线 可 以 有 效 地 反应 不 同 区 域内 气象 条 件 以 及 
自然 环境 的 差异 。 根据 在 Bagrot 山谷 观测 期 间 所 获 
得 的 降水 样品 的 8D F SO 数据 ,利用 最 小 二 乘法 ， 

分 别 计算 得 出 研究 区 全 年 局 地 大 气 水 线 (图 4) : 


8D =7.778"0O -4.28(R^ =0.99 =83,P<0.01) 
春 夏季 大 气 水 线 : 

8D 27.8980 +3. 82(R° 20.98,n 258,P «0.01) 
秋冬 季 大 气 水 线 : 

8D 27.8180 +6.80 (R^ 20.99,n 225,P «0.01) 


与 全 球 大 气 水 线 相 比 , 全 年 局 地 大 气 水 线 斜率 
和 截 距 均 偏 低 ,这 可 能 是 由 于 降水 过 程 中 雨滴 受到 
了 较为 强烈 的 云 下 二 次 落 发 的 作用 引起 的 ,也 即 间 
接地 受 控 于 当地 相对 干燥 的 自然 环境 。 已 有 研究 表 
8] , 当 气 温 较 高 且 相 对 湿度 较 低 的 情况 下 ,降雨 过 程 
所 受到 的 云 下 二 次 蒸发 作用 较 强 ,由 此 造成 了 同位 
RAE IME OE 。 研 究 时 段 内 较 小 的 降水 事 
件 (<5 mm) 占 总 降水 事件 次 数 的 65% ,降水 量 较 
小 的 降水 事件 更 容易 受到 云 下 二 次 蒸发 的 影响 ,从 
而 导致 研究 区 大 气 水 线 斜率 和 截 距 低 于 全 球 大 气 水 
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图 4 Bagoi 山谷 全 年 ,不 同 季节 ( 春 夏季 、 秋 冬季 ) 的 
局 地 大 气 水 线 和 全 球 大 气 水 线 
Fig.4 GMWL and LWML for different seasons 


( spring-summer , autumn-winter) in Bagrot Valley 


线 , 这 一 过 程 在 较为 干旱 的 地 区 尤为 强烈 ” S 53 
外 ,研究 还 发 现 , 降 水 的 不 同形 态 对 局 地 大 气 水 线 的 
斜率 和 截 距 产生 一 定 的 影响 。 本 研究 的 所 有 降雨 事 
fF 8D 和 5"0 之 间 的 线性 关系 (图 5) 为 : 


8D 27.8280 43.99 (R? 20.98,n 271,P «0.01) 


相对 而 言 ,降雪 过 程 则 会 造成 截 距 和 和 斜率 的 偏 
大 ,本 研究 中 所 有 降水 样品 中 的 12 个 降雪 样品 ,其 
8D fi 8^0 之 间 的 线性 关系 (图 5) 为 : 


8D 29.0080 +26. 27( R^ =0.99,n=12,P<0.01) 


式 中 的 斜率 和 截 距 均 大 于 全 球 大 气 水 线 的 以 及 
研究 区 全 年 大 气 水 线 的 ,表明 冬季 降雪 事件 受 再 循 
环 水 汽 和 二 次 蒸发 的 影响 较 小 ,水 汽 来 源 较为 单一 。 
这 也 是 秋冬 季 局 地 大 气 水 线 比 春 夏季 局 地 大 气 水 线 
和 斜率 和 截 距 大 的 原因 。 

而 通过 与 研究 区 周边 地 区 局 地 大 气 水 线 对 比 发 
现 ( 表 3) ,整体 上 研究 区 周边 地 区 大 气 水 线 与 研究 
区 大 气 水 线 相 比 ,其 斜率 大 小 比较 相近 ,但 其 截 距 显 
著 不 同 。 研 究 发 现 , 从 南 到 北 , 总 体 上 随 着 纬度 的 升 
we 
除外 ) 。 这 可 能 是 由 于 随 着 纬度 的 升 高 ,降水 受 海 
影响 越 小 ,气候 愈加 干燥 雨滴 降落 时 强 
烈 的 云 下 二 次 莹 发 作用 使 得 稳定 同位 素 发 生 不 平衡 
分 馏 , 造 成 截 距 逐 渐 降 低 。 

2.3 水汽 来 源 分 析 

从 地 理 位 置 上 来 看 ,研究 区 远离 海洋 ,位 于 亚洲 

大 陆 内 部 ,并 且 研 究 区 南部 有 众多 山脉 ,阻挡 了 大 部 
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图 5  Bagrot 山谷 降雨 .降雪 不 同 降水 形态 时 的 
局 地 大 气 水 线 与 全 球 大 气 水 线 


Fig.5 GMWL and LWML for different precipitation phases 


(rainfall events and snow events) in Bagrot Valley 
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RI 研究 区 以 及 周边 地 区 局 地 大 气 水 线 ( 以 纬度 排序 ) 


Tab.3 Local Water Meteoric Line of Bagrot Valley and its adjacent area 

研究 区 域 经 度 纬度 局 地 大 气 水 线 资料 来 源 
ANE 79. 58° 37. 10° òD =7.248'%0 +1. 96 GNIP 
Bagrot Valley 74. 56° 36. 03^ 8D 27.7780 +4. 28 本 研究 
Kashmir Himalaya 74. 50° 34.25° 8D =7.598'80 +11.79 文献 [49 -50] 
阿里 79.70° 33. 39° 8D =8.518'°0+11.57 文献 [56 ] 
狮 泉 河 84.08。 32. 50° 8D -8.098/50 +7.96 文献 [42 ] 
改 则 84. 02° 32. 03° 8D 28.2980 +11.62 文献 [57 ] 


分 来 自 于 印度 详 的 暖 湿 气 流 ,降水 可 能 主要 来 自 于 
西风 带 气流 的 水 汽 输 送 以 及 局 地 环流 。 根 据 YAO 
等 人 的 研究 结果 ,可 以 认为 本 研究 区 域 降水 隶属 于 
西风 模 态 ,降水 稳定 同位 素 变 化 主要 受 温度 变化 控 
W, 

为 进一步 研究 该 区 域 水 汽 来 源 ,借助 NOAA 提 
供 的 HYSPLIT 气 团 轨迹 模型 ,对 Bagrot 山谷 样品 收 
集 站 点 上 空 1 500 m 处 的 气 团 为 运动 终点 ,对 2015 
年 8 月 到 2016 年 7 月 降水 进行 后 向 120 h 轨迹 模拟 
追踪 ,将 模拟 轨迹 近似 代表 水 汽 团 的 轨迹 "| 。 

春 夏 季 和 秋冬 季 降 水 的 水 汽 运行 轨迹 有 一 定 差 
别 ( 图 6a,6b) ,可 以 将 春 夏 和 秋冬 季 的 水 汽 来 源 大 
致 分 为 3 大 类 :(1) 来 自 于 研究 区 西部 的 大 西洋 、 亚 
欧 大 陆 腹 地 以 及 地 中 海区 域 ;(2) 来 自 于 研究 区 周 
边 地 区 ;(3) 来 自 于 研究 区 南部 的 阿拉 伯 海 。 

第 一 类 水 汽 来 源 主要 经 西风 环流 长 距离 输送 到 
达 研 究 区 ,地 中 海地 区 和 亚 欧 大 陆 腹 地 为 降水 d 的 
峰值 区 "1 , 故 来 自 于 此 区 域 的 水 汽 具 有 和 较 高 的 d 
值 ,可 以 看 到 秋冬 季 由 西风 环流 来 的 水 汽 轨迹 距离 
相 较 于 春 夏季 的 更 长 ,并 且 占 总 轨迹 数 的 比例 更 高 ， 
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注 : 秋 冬季 为 9 月 ~2 月 ; 春 夏季 3 月 ~8 月 ;轨迹 计算 为 在 降水 发 9 


导致 4 值 偏 高 而 3 "0 值 偏 低 。 第 二 类 水 汽 来 源 为 
局 地 循环 水 汽 , 通 过 短 距离 的 水 汽 传输 过 程 到 达 研 
究 区 ,此 类 水 汽 一 般 来 自 于 土壤 蒸发 .植物 蒸腾 以 及 
地 表 水 体 蒸发 , 春 夏季 由 于 温度 较 高 ,蒸发 较 强 ,对 
流 天 气 较 多 ,强烈 的 再 循环 过 程 使 得 5^0 富 集 , 并 
fi d 值 偏 低 '”] 。 第 三 类 水 汽 来 源 主要 为 来 自 南方 
的 海洋 性 水 汽 。 此 类 水 汽 从 沿海 到 内 陆 经 历 长 途 输 
送 ,在 输送 过 程 中 水 汽 不 断 “ 洗 脱 ”" 以 及 遇 到 高 大 山 
脉 阻 挡 被 迫 抬 升 , 从 而 凝结 产生 降水 ,使 得 剩余 水 汽 
8^0 不 断 贫 化 ,到 达 研 究 区 后 ,导致 后 续 形成 的 降 
水 出 现 较 低 的 3 O 值 ,这 可 以 解释 研究 区 偶尔 出 现 
的 较 低 的 8*0 的 降水 过 程 ”1。 例 如 3 月 18 日 ,由 
HYSPLIT 模型 的 模拟 结果 来 看 (图 7) ,其 后 向 轨迹 
追踪 在 3 000 m 高 度 有 来 自 于 阿拉 伯 海 的 水 汽 输 
入 ,从 而 导致 Bagrot 山谷 此 日 降水 稳定 同位 素 相 对 
较 低 ,3850 和 8D 分 别 为 - 13. 83%o 和 — 113. 1596. 
但 3 月 19 日 模拟 结果 显示 ,水 汽 主要 来 自 于 西风 远 
距离 和 局 地 循环 ,导致 此 日 降水 稳定 同位 素 相对 较 
$3,860 和 8D 分 别 为 -10.03%o 和 一 70.57%o。 可 见 
不 同 水 汽 来 源 的 转换 会 影响 降水 稳定 同位 素 的 变化 。 
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事件 当天 的 0:00 ,6 00 ,12 00 Fil 18 :00 ,模拟 时 间 为 向 前 120 h, 


高 度 为 地 面 以 上 1 500 m ,数据 为 NCEP 再 分 析 数 据 
图 6 Bagrot 山谷 气流 后 向 轨迹 图 
Fig.6 Backward trajectories simulated by HYSPLIT for Bagrot Valley 
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注 :红色 轨迹 表示 地 面 上 模拟 高 度 1 000 m; 蓝 色 轨 迹 表 示 地 面 上 2 000 m; 绿 色 轨 迹 表示 地 面 上 3 000 m 


7 HYSPLIT 后 向 轨迹 追踪 模拟 图 
Fig.7 Back trajectories simulated by HYSPLIT at 1000 (red lines) , 2 000 (blue line) and 3 000 m (green line) a. g.l. on 
March 18 and 19 ,2016 


从 整体 上 看 ,Bagrot 山谷 降水 稳定 同位 素 变化 
村 征 与 青藏 高 原 西 北部 塔 什 床 尔 干 的 研究 结果 相 
似 , 其 变化 主要 受 西风 环流 和 局 地 循环 水 汽 的 影 
jj^" 。 不 过 ,与 塔什库尔干 以 及 和 田 等 地 的 研究 稍 
有 不 同 的 是 ,Bagrot 山谷 降水 稳定 同位 素 仍然 识别 
到 南方 的 海 详 性 水 汽 偶尔 仍 会 影响 到 此 研究 区 。 这 
可 能 与 Bagrot 山谷 相 比 塔什库尔干 更 靠 南 ,并 且 海 
拔 较 低 有 关 。 


3 结论 


本 研究 通过 对 印度 河上 游 流 域 Bagrot 山谷 
2015 年 8 月 一 2016 年 7 月 降水 稳定 同位 素 观测 结 
Ae ,结合 研究 区 气象 观测 数据 以 及 HYSPLIT 模型 后 
向 轨迹 追踪 ,利用 同位 素 示 踪 及 统计 分 析 方 法 分 析 
了 Bagrot 山谷 降水 稳定 同位 素 变化 特征 及 与 水 汽 来 
源 的 关系 ,得 出 以 下 主要 结论 : 

(1) Bagrot 山谷 降水 稳定 同位 素 8^0 5j se 
节 变 化 明显 ,秋冬 季 较 低 , 春 夏季 较 高 ,与 温度 呈正 
相关 关系 ,全 年 温度 效应 明显 ,相关 系数 达到 了 
0.73 ,而 降水 量 效应 不 明显 。 

(2) Bagrot 山谷 局 地 大 气 水 线 截 距 和 和 斜率 均 小 
于 与 全 球 大 气 水 线 , 这 反映 了 研究 区 较为 干旱 的 环 
境 和 降水 受 云 下 二 次 蒸发 的 影响 。 特 别 的 是 ,通过 
对 不 同 的 降水 形态 (降雨 .降雪 ) 大气 水 线 研 究 发 


现 , 不 同 的 降水 形态 会 导致 该 研究 区 的 局 地 大 气 水 
线 的 斜率 和 截 中 不 同 。 降 雨 事件 受到 的 二 次 燕 发 相 
对 较 强 ,使 得 局 地 大 气 水 线 的 斜率 和 截 距 明显 偏 低 ; 
而 降雪 事件 发 生 时 ,由 于 环境 温度 较 低 , 受 再 循环 水 
汽 和 二 次 蒸发 的 影响 较 小 ,导致 其 局 地 大 气 水 线 的 
斜率 和 截 距 偏 高 。 通 过 与 研究 区 周边 地 区 局 地 大 气 
水 线 对 比 发 现 ,从 南 到 北 ,其 斜率 大 小 基本 相同 ,而 
其 截 距 随 着 纬度 的 升 高 呈 变 低 的 趋势 ( t TRES a 
除外 ) ,这 可 能 是 由 于 随 纬度 升 高 , 受 海洋 性 水 汽 影 
响 越 小 , 受 云 下 二 次 藻 发 影响 越 大 造成 的 。 

(3) 根据 降水 稳定 同位 素数 据 及 HYSPLIT 模 
型 模拟 结果 表明 :研究 区 全 年 主要 受 西风 环流 以 及 
局 地 环流 的 影响 ,不 同 水 汽 来 源 的 转换 会 引起 降水 
稳定 同位 素 的 变化 。 但 是 ,与 研究 区 以 北 的 临近 站 
点 ( 荣 土 塔 格 、 和 田 等 ) 相 比 有 所 不 同 ,由 于 Bagrot 
山谷 位 置 更 靠 南 ,其 仍然 偶尔 受到 来 自 南 方 的 海洋 
性 水 汽 影响 。 
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Abstract: Precipitation stable isotopes are integrated tracers of atmospheric processes worldwide. The precipitation 
stable isotopes (80 and 8D) are influenced by numerous factors , including the condensation temperature , atmos- 
pheric humidity , air pressure, precipitation amounts, moisture sources and orographic terrain. In recent decades, 
many efforts have been dedicated to studying precipitation isotopic composition in various regions in the Tibetan 
Plateau , China, located at the convergence between the westerlies and Indian monsoon. However, few studies are 
conducted in the Upper Indus Basin ( UIB) , which located in the western Tibetan Plateau. Based on the observation 
of precipitation stable isotopes (8^0 and 8D) and local meteorological factors from August 2015 to July 2016 , com- 
bining with the HYSPLIT model ,this study investigated the characteristics of stable isotopes in precipitation ,the re- 
lationship between isotopic variations and meteorological factors, Local Meteoric Water Line ( LMWL) , and the 
moisture sources of typical watershed in Bagrot Valley , Upper Indus Basin. The results show that precipitation 8^0 
and 8D in Bagrot Valley displayed an obvious seasonal change. The 8^ O and 8D values are more depleted in au- 
tumn and winter but relatively enriched in spring and summer , indicating a significant temperature effec rather than 
precipitation amount effect. Moreover ,the lower intercept and slope of Local Meteoric Water Line (LMWL) relative 
to GMWL indicated a strong sub-cloud evaporation owing to the relatively arid climate in the study area. Interesting- 
ly , different precipitation phases led to a difference in intercept and slope of meteoric water line. The lower intercept 
and slope of rainfall events result from strong sub-cloud evaporation. While meteoric water line of snowfall events 
has higher slope and intercept. It resulted from little impacts of recycled moisture and sub-cloud evaporation , due to 
the low temperature. Besides , compared with adjacent regions ,the slopes of the LMWLs are nearly the same , howev- 
er, the intercepts were increased with the latitude (except Shiquanhe) ,this may be caused by the strengthening sub- 
cloud evaporation and weakening effect of the marine moisture. In addition , the outputs of HYSPLIT model and re- 
sults of precipitation stable isotopes demonstrated that Bagrot Valley is consistently dominated by the westerlies and 
local moisture recycling throughout the year and moisture source changes affect stable isotopes in precipitation. How- 
ever,compared with Muztagata and Hotan , where marine moisture cannot reach, the study area of Bagrot Valley is 
affected by the marine moisture from the ocean occasionally as it lies southern to Muztagata and Hotan. 
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